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ABSTRACT 

Gravitational wave emission is considered to be the driving force for the evolution of short- 
period cataclysmic binary stars, making them a potential test for the validity of General Rel- 
ativity. In spite of continuous refinements of the physical description, a 10% mismatch ex- 
ists between the theoretical minimum period (Ptum — 70 min) and the short-period cut-off 
(P m ; n ~ 80 min) observed in the period distribution for cataclysmic variable binaries. A pos- 
sible explanation for this mismatch was associated with the use of the Roche model. We here 
present a systematic comparison between self-consistent, numerically constructed sequences 
of hydrostatic models of binary stars and Roche models of semi-detached binaries. On the 
basis of our approach, we also derive a value for the minimum period of cataclysmic variable 
binaries. The results obtained through the comparison indicate that the Roche model is indeed 
very good, with deviations from the numerical solution which are of a few percent at most. 
Our results therefore suggest that additional sources of angular momentum loss or alternative 
explanations need to be considered in order to justify the mismatch. 
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1 INTRODUCTION 

Semi-detached binary systems composed of a white dwarf as the 
primary and a low-mass main sequence star or brown dwarf as 
the se condary, are u sually referred to as Cataclysmic Variables, or 
CVs (Warner 1995J). Such systems could show a cataclysmic be- 



haviour as a consequence of mass overflow from the secondary to 
the primary star. The observed orbital period distribution of these 
systems has two distinctive features: (a) a statistically significant 
deficit of systems with orbital periods between 2 and 3 hours, 
known as the "period gap"; (b) a sharp lower cutoff at an orbital 
period o f -Pmin — 80 min w hich is referred to as the "minimum 
period" (JRitter & Kolb 1998h. 



In modeling short-period CVs (i.e. CVs with an orbital pe- 
riod < 2 hr), it has long been proposed that gravitational radiation 
reaction drives their evolution enabling a continuous mass flow to 
occur from the secondary to the primary star (Kraft, Mathews & 
Greenstein 1962; IPaczyriski 198 1|: |Paczyhski & Sienkiewicz 1981 ; 



Rappaport, loss & Webbink 1982|; P'Antona & Mazzitelli 1982 



Paczyhski & Sienkiewicz 1983|; Nelson, Chau & Rosenblum 1985, 



Kolb & Ritter 1992|; |Howell, Rappaport & Politano 1997J Kolb~& 
Baraffe 19 99). This feature makes CVs a p otential test of General 
Relativity ( Paczynski & Sienkiewicz 1983| 



Note that for CVs with an orbital period > 3 hr the emission of gravita- 
tional waves cannot drive the binary evolution. 



The evolutionary scenario for a short-period CV can be sum- 
marized as follows. Consider a binary system with the secondary 
star having a mass O.1M < M 2 < 0.3M Q Q As a result of the 
binary evolution, the system becomes semi-detached, with the sec- 
ondary star filling its "critical lobe" (which is defined here as the 
three-dimensional equipotential surface having a cusp at the inner 
Lagrangian point in The loss of angular momentum then induces 
a mass overflow from the inner Lagrange point onto the primary 
degenerate dwarf. For the particular class of CVs known as Dwarf 
Novae an accretion disc is formed around the primary dwarf star 
and it is this which is responsible for the observed "cataclysmic" 
behaviour. 

During the initial stages of the binary evolution, the mass- 
transfer timescale r^, which is of the same order as the timescale 
for angular momentum loss due to the gravitational wave emission 
tgr, is much longer than the Kelvin-Helmholtz timescale r KH . Un- 
der these conditions, the secondary can maintain its thermal equi- 



T A slightly different evo lution scenario should be considered in the case 
in which M 2 < O.IMq ( paczynski 1981 ). 

■f We here introduce the term "critical lobe" in order to distinguish be- 
tween the critical equipotential surface obtained from a self-consistent 
three-dimensional solution of the hydrodynamical equations and the corre- 
sponding surface obtained when using the Roche model and which is usu- 
ally referred to as the "Roche-lobe". 
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librium and it slowly moves down the main sequence as it loses 
mass. As a result of angular momentum loss, the orbital period 
of the binary system decreases and the mean density of the sec- 
ondary increases. A direct consequence of this is that r KH for the 
secondary also increases. 

Because in the course of the evolution the ratio Tj^/tkh de- 
creases secularly, the secondary gets more and more out of thermal 
equilibrium. This is connected with a tendency of the star to shrink 
less and less or even to expand upon mass loss, depending on the 
value of Tj^/tkh. The minimum (turning) period of the binary or- 
bit Ptum is reached for a specific value of the ratio Tj^/tkh which 
depends somewhat on the model assumed for the secondary [cf. eq. 
(2)]. After the minimum period is reached, the nuclear burning be- 
comes inefficient and radiative losses from the surface of the star 
will reduce its heat content. Eventually the star expands despite 
losing mass and its evolutionary path moves away from the main 
sequence towards the degenerate star branch. 

Within this evolutionary scenario, Ptum is usually determined 
by means of the so called Roche model, which we will discuss in 
more detail in Section 3.1 . The minimum period found b y several 
authors in this way (Paczyhski 1981; JKolb & Ritter 1992[ Howell, 



Rappaport & Politano 1997), turns out to be about 10 % shorter 
than the observed period cut-off P m i n (see Kolb & Baraffe 1999 
for the latest theoretical results). This difference seems to survive 
all attempts at introducing refined treatments of the stellar physics 
and of the stellar evolution. 

There have been many proposals for solving this discrepancy 
between Ptum and P m in. One of the suggestions is that the finite 
size effects of the secondary star, that are neglected in the tradi- 
tional Roche model, ma y not be negligible and might be re spon- 
sible for the m ismatch (Nelson, Chau & Rosenblum 1985[ K olb 
& Baraffe 1999; |Lasota 2000| ). Indeed, Nelson et al. (1985) have 
claimed that when tidal distortion effects are included in their mod- 
els on the basis of the scheme suggested by Chan & Chau (1979), 
Ptum is increased by ~ 10 %. These fi ndings, however, were not 
confirmed by subsequent investigations ( ^Colb & Baraffe 1999 ). 

In this paper, we focus our attention on the validity of the 
Roche model for calculating Ptum in CVs. For this doing, we com- 
pare the dynamical properties of semi-detached binaries obtained 
by means of the Roche model with those obtained by numerically 
solving the equations of hydrostatic equilibrium for a binary sys- 
tem. In our numerical approach, we solve exact equations for the 
hydrostatic equilibrium, so that the tidally deformed structure of 
the secondary star and the orbital motion of the binary system are 
computed self-consistently. This is a feature missing in the Roche 
model. In order to emphasize the differences which emerge from 
the hydrodynamics of the equilibrium binary system, our calcula- 
tions do not include details of the thermal structure or evolution of 
the secondary. 

The paper is organized as follows: in Section we describe the 
assumptions made in the construction of the models for the binary 
system. In Section pi we discuss the set of equations that is used in 
the numerical construction of the equilibrium sequences of binary 
stars. Sections Wand ra summarize the numerical results and present 
our conclusions. 



2 ASSUMPTIONS OF THE MODELS 

A number of assumptions are necessary for the construction 
of a model for a CV binary system. In particular, following 
Paczyhski (1981), we consider a semi-detached binary system with 



a white dwarf as the primary star with mass Mi, and a hydro- 
gen rich star on the lower main sequence as the secondary star 
with mass Mi- Mi and M2 are assumed to be typically Mi ~ 
0.5 - 1.0M© and M 2 < O.4M . Since the secondary is taken to 
be fully convective in this mass range, its structure is well described 
by a polytropic relation between the pressure p and the mass den- 
sity p 



p = K P 1+ t? 



(1) 



with polytropic index N = 3/2 and polytropic coefficient K, 
which is constant in space for an isentropic modelpl For such a 
polytrope, the condition for the occurence of the minimum p eriod 
can the n be expressed as (fr)' Antona, Mazzitelli & Ritter 1989; Rit- 
ter 1996; [Kolb & Ritter 1992J ) ~ 

^- ~ 0.3 . (2) 

This condition, which contains some uncertainties related to the 



stellar model assumed, will be used in Section 12 to determine the 



range of minimum periods from the numerically calculated values 
of the mass transfer timescale t^. 

The radius of the primary is small compared to that of the sec- 
ondary and the primary is therefore treated as a point-like mass. 
This simplifies the calculation of both the gravitational potential of 
the primary and its contribution to the quadrupole moment of the 
system [cf. eq (|2l|)]. The binary orbit is assumed to be circular and 
the spin period of the secondary is taken to be synchronous with 
the orbital period. These are rather good approximations since the 
relaxation timescales to these state s are much shorter than t he evo- 
lutionary timescale of t he system (|Zahn 1975J |Zahn 1977|; L ecar, 
Wheeler & McKee 1976). 

Finally, we assume that the secondary is always just filling its 
critical lobe and that the continuous mass-transfer from the sec- 
ondary to the primary is conservative, so that the total mass of the 
system M = Mi + M2 is constant in time. This latter assumption 
is justified when we consider the evolution of the binary system to 
be driven solely by the emission of gravitational radiation. 



3 CONSTRUCTION OF THE MODELS 

In this section we briefly discuss the two distinct sets of equations 



3.1 



for the construction of models for the CVs. In Section '. 

marize the set of equations describing the "standard" Roche model 

and in Section 3.2, we discuss the set of the equations for a numer- 



ically constructed hydrostatic model of a semi-detached binary. 



3.1 Roche model of a semi-detached binary 

We construct Roche models of CVs based on the Roche model of 
semi-detached binaries as discussed by Paczyhski in 1981. In this 
approximation, the secondary star is assumed to fill up its Roche 
lobe and its gravitational field is treated as that of a point-like mass. 
Within this approach, it is convenient to measure the size of the 
Roche lobe through a mean radius r L defined so that the sphere of 
radius r, has the same volume Vo as the Roche lobe 



(£)' 



(3) 



3 Note that K changes with q along the sequence of constant total mass 
equilibrium models. 



© 0000 RAS, MNRAS 000, 000-000 



Gravitational Wave Emission in CVs 



The ratio of r L to the binary separation a depends only on the 
mass ratio q = M2/M1 and the empirical expressions proposed by 
Paczyhski (1971) and Eggleton (1983) are expressed as 



f(l) 



2 

it 



0.49g3 



0.6g? +ln(l + q*) 



(Paczyhski 1971) 



(Eggleton 1983) 



(4) 



The above expressions are valid for the range of mass ratios of in- 
terest here: q ~ 0.1 — 1, Paczyhski's expression is accurate to 
within 1% when compared wit h the result o btained by numerical 
integration of the Roche-lobe (Kopall959). In the same range, 
Eggleton's expression is even more accurate. 

Combining the above formula with Kepler's law 



a = (£) G(M 1 + M 2 ) , 

we obtain an expression for the orbital period P as 



P = 



2tt 



M 2 



(l + q)(f(q)) 3 



(5) 



(6) 



Within the Roche model, the spin angular momentum of the pri- 
mary star is neglected and the total angular momentum J of the 
system is therefore given by 



'J — 'Jorh T </spin 



(7) 



where the orbital angular momentum of the system J or b and the 
spin angular momentum of the secondary J sp i n are defined respec- 
tively as, 



•-'or 



and 

'/soil 



Ga 



' M ± + M 2 



= G? 



m 



q(l + q) 



(r g r L yM 2 n = (1 + q)f(q)r 2 g J OTh . 



(8) 



(9) 



Here r g rL is the radius of gyration of the secondary and for a poly- 



trop e with N = 3/2, 



_ 



0.45 (D'Antona, Mazzitelli & Ritter 



1989). The theory of stellar structure provides a simple mass-radius 
relation which holds for lower main sequence stars in the mass 
range of interest here (see, for example, Kippenhahn & Weigert 
1990). This relation can be written as 



R VJ 



X. 



M 2 
M Q 



(10) 



where Mq and Rq are the solar mass and the solar radius, respec- 
tively, and R 2 is the radius of the lower main sequence star. For a 
star with chemical abundances similar to those of the Sun, we can 
set x = 0.83 and £ = 0.84 (Kippenhahn & Weigert 1990). More- 
over, as long as the star remains on the lower main sequence (i.e. as 
l ong as r KH < r^), we can assume that the mass-radius relation 
(jlOp holds even in the presence of mass loss, i.e. when 



R 2 



(11) 



Using equations (|10J) and (|ll|), we can now rewrite the expression 
for the orbital period (n) as 



P = 2.78 x 



35-1 



M-> 



Mq (l + g)(/(g))3 



hr , (12) 



and derive the two timescales which are relevant for the evolution 
of a CV binary: namely the timescale for gravitational wave emis- 
sion tgr and that for mass-transfer r^. As mentioned above, we 
assume that the binary evolution is driven only by the emission of 
gravitational waves (GW) and the timescale for this to occur can 
be estimated by means of the standard quadrupole formula for two 
poin t-like masses orbiting around each other (L andau & Lifshitz 
1975) - 

din J ^ - 1 



/ d In J \ 

V dt gw/ 

32 \2irJ 



3 c 5 (Mi + M 2 ) 
G § M1M2 



-7.87 x 10 



'(1 



hi J 



M 2 

M 



[q^l + q)] 1 yr, 



(13) 



On the other hand, we can calculate the mass-transfer timescale as 

'dlnM 2 \- 1 fdlnAfa din J ■> _1 



'M 



lit 



GW 



(14) 



d In J dt 

where the logarithmic differential of the total angular momentum 
am can be expressed as 



din J = 



1 3 

-dlnr L +-dlnM 2 



1 



din Mo 



1 

~2 


'din /(g) 2g + l " 
dg g(l + g)_ 


dq 


2-2g-g(l + g)^ln/(g)+£ 



,(15) 



where we have used equations (JlOJ) and (hill as well as the conser- 
vation of the total mass of the system 



dg = g(g+l)dlnM 2 . (16) 

Combining equations ( |13p - (Jl5[), we can write explicitly the mass- 
transfer timescale for the secondary as 

dlnJ X- 1 

GW/ 



M\A 



l( : 



(It 



2q - q(l + q)^- In f(q)+C 



(17) 



Note that the CV system described by the Roche model is deter- 
mined in terms of two variables: the mass ratio q and the mass of 
the secondary, M 2 . 

3.2 Numerical hydrostatic model of a semi-detached binary 

A numerical method for computing semi-detached compressible bi- 
naries in hydrostatic equilibrium has recently been developed by 
Uryu & Eriguchi (1999) [see also Uryu & Eriguchi (1998) for 
the case of incompressible binaries.]. Within the assumptions dis- 
cussed in Section Q, and in a corotating reference frame, the matter 
of the secondary star satisfies the following equation of hydrostatic 
equilibrium (the Bernoulli equation) 



-zu 2 0. 2 + K(N + 1)9 + </> = C . 



(18) 



where, <p and Q. are the gravitational potential, and the angular ve- 
locity of the orbital motion, respectively, w is the distance of the 
fluid element from the rotation axis of the binary and C is a con- 
stant. The function S is defined by p = . The gravitational 
potential <f>, is computed using the integral expression 
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-G 



P(r') 



dV 



Ghh 



(19) 



V-2 



where the volume integration is carried out in the interior of the sec- 
ondary (V2), and ri is the position vector of the primary. Equations 
< jl 8[ > and the integral in the equation ( |19| ) are discretized in a spher- 
ical coordinate system (r, 9, ip) whose origin is at the centre of the 
secondary. Because of the symmetry of the configuration, we only 
compute the solution in the region r G [0, Ro], 9 G [0, 7r/2], tp G 
[0, 7r], where Ro is the maximum radius of a tidally deformed stel- 
lar surface. A solution is obtained on a grid with the number of 
gridpoints for each coordinate being (N r , Ng,N v ) — (65, 33, 65). 
The Green's function of the integral in the equation (fl9) is ex- 
panded in Legendre polynomials Pi (cos 9) up to £ — 14. The num- 
ber of grid points as well as the number of Legendre polynomials 
used for the Green's function are chosen so as to provide satisfac- 
tory numerical accuracy and results presented in this paper have a 
numerical error which is typically less than 0.1%. [Further details 
of the numerical method can be found in Uryu & Eriguchi (1999) 
and references therein.] 

In calculating models of "critical lobe-filling binaries" (here- 
after critical lobe binaries), we first fix a value for the mass ratio q 
and then compute a set of solutions corresponding to different val- 
ues of the orbital separation a (we start from a large initial value of 
a and then decrease it in steps of 8a — Ro/50). The sequence of 
solutions obtained in this way is terminated at an orbital separation 
for which the surface of the secondary star has a cusp at the in- 
ner Lagrange point. The corresponding configuration is identified 
as a critical lobe model for the given value of q. This procedure 
has been repeated for 3 1 models, changing q logarithmically in the 
range gG [0.033, 1]. 

For each of the 31 models, we compute the relevant timescales 
and, as in equation dob, r GW is given by the rate of angular mo- 
ment um loss from the system in the form of gravitational radia- 
tion (Misner, Thorne & Wheeler 1975 ), 



din J 



dt 



GW 



32G» 5 
""5c5"T 



(/11-/22) 2 



(20) 



The quadrupole moment tensor hj of the system is computed nu- 
merically as 



(Ps +Pp) 



dV 



(21) 



where the volume V comprises both stars and p s , p p are the mass 
densities of the secondary and primary stars, respectively, with the 
latter being expressed in terms of a three-dimensional Dirac delta 
function p p = Mi8 a (y — zuie xi ). The vector r = (x\,X2,%z) 
is the position vector of a mass element and wi represents the dis- 
tance of the primary star (a point particle) from the rotation centre 
of the binary. The vector e xi is the unit vector along the x\ axis. 
The Latin indices in (Ell) run from 1 to 3 and correspond to the com- 
ponents of the Cartesian coordinate system whose origin is at the 
rotation centre of the binary system (the x%— coordinate is oriented 
along the separation vector of the binary and the X2— coordinate is 
in the orbital plane). The total angular momentum J for the model 
is computed as follows : 



J 



pvu Q,dV + Miuj 1 0, . 



(22) 



V"2 



As in equation (M), we calculate the mass-transfer timescale as 
'dlnM 2 \' 1 /dlnAf 2 din J 



m -c 



d In J dt 



(23) 
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Figure 1, The relative differences of several quantities when computed us- 
ing the Roche model and the numerical binary models are plotted against 
q. The solid line is the result for r, as computed with Paczyhski's formula. 
The short-dashed line is the same but for Eggleton's formula. The long- 
dashed line, is the relative error for (In — I22) 2 [see eq. (H)]- 



where now the quantity din M2/ din J is computed along the 
semi-detached equilibrium sequence from the expression 

dlnM 2 1 dq 



din J 



q(q+ 1) din J 



(24) 



4 COMPARISON OF THE MODELS 

In this Section, we compare the two classes of semi-detached bi- 
nary models, namely the Roche model and the numerical critical 
lobe model in hydrostatic equilibrium. We first compare purely dy- 
namical quantities for the two models. We then discuss the differ- 
ences arising from the application of these two models to a cata- 
clysmic binary system, focusing in particular on the differences in 
the mass-transfer timescale. 



4.1 Comparison of dynamical properties 

Two binaries are considered to represent the same physical system 
(and are therefore considered "comparable") if they have the same 
mass ratio q and the same orbital period P. In practice, we first 
construct a numerical model of a critical lobe binary and then cal- 
culate the physical quantities of Roche model binary with the same 
mass ratio and orbital period. Note that the relation (flfJ) between 
the mean radius and the mass of the secondary is not used in this 
comparison, but we match two different binaries on the basis of 
the dynamical timescale of the system (i.e. the orbital period of the 
binary system). 

In the comparison we focus on: (a) the mean radius of the 
secondary r L ; (b) the total angular momentum of the system J; 
(c) the components of the quadrupole moment tensor (In — I22) 2 
and (d) the inverse of the timescale for gravitational wave emission 
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Figure 2. The relative deviations of the total angular momentum when com- 
puted using the Roche model and the numerical binary models are plotted 
against q. The solid line refers to binaries in which the spin angular momen- 
tum of the secondary is neglected in the Roche model. Binaries in which the 
spin contribution of secondary in the Roche model is instead included are 
indicated with a short-dashed line (for the case of Paczyhski's formula) and 
with a long-dashed line (for the case of Eggleton's formula). 



Figure 3. Relative error of the inverse of the timescale for angular mo- 
mentum loss by gravitational radiation, d In J/dt. The solid line refers to a 
binary system with total mass M = 0.8Mq and does not include the spin 
angular momentum contribution of the secondary, which is instead included 
in the curve indicated with a short-dashed line. Similarly, the long-dashed 
and dotted lines are for a system having a total mass M = 1.3Mq, with 
and without the inclusion of the spin angular momentum of the secondary. 



d In J/dt. For all of these quantities we define the measure of the 
relative deviation as 



Relative error in Q = 



Q (critical lobe) — Q (Roche model) 
Q (Roche model) 



(25) 

and have plotted them in Figures [lj-pj against the mass ratio of the 
system, and in Figure against the orbital period. 

In particular, the solid line in Figure hi shows the mean radius 
of the critical lobe, r L , from Paczyhski's formula, while the short- 
dashed line comes from using Eggleton's formula. The long-dashed 
line, on the other hand, shows the relative error for (In — /22 ) 2 P| 
Similarly, Figured shows the relative deviations of the total angular 
momentum: the solid line refers to binaries in which the spin angu- 
lar momentum of secondary is neglected in the Roche model, while 
binaries for which the spin contribution of secondary in Roche 
model is included are indicated with short and long-dashed lines. 
Both Figures m and PJ show that for all of these quantities the errors 
made through the use of the Roche model are smaller than 3%, if 
the spin angular momentum of the secondary is taken into account. 

In a similar way, in Figure H, where we calculate the rela- 
tive error in d In J/dt for different values of the orbital period, 
the different line types refer to two systems having a total mass 
of M = 0.8M Q and M — 1.3M , respectively. In this compari- 
son we have used Eggleton's formula but we have verified that the 
error is almost the same when Paczyhski's formula is used. Also 
in this case, the deviations are of at most a few percent, as found 
by Kolb & Baraffe (1999) when implementing the prescription of 
Chan & Chau (1979) for the corrections due to tidal distortion. 

^ We recall that in the Roche model h 1-/22 = a 2 M\Mi / (M\ + A-/2). 
© 0000 RAS, MNRAS 000, 000-000 



Note that in evaluating the timescale for the angular momen- 
tum loss, the accuracy of the Roche model increases as the mass 
of the secondary is decreased and thus for shorter orbital periods. 
This can be deduced from the fact that, for a fixed value of the total 
mass M, the relative error in d In J/dt becomes smaller for smaller 
orbital periods and therefore for smaller masses of the secondary. 
The improved accuracy is related to the fact that the Roche model 
is increasingly accurate as the finite-size effects become less and 
less important, or equivalently as the ratio r L /a becomes smaller 
and smaller. In the case of Paczyhski's formula, this ratio can be 
expressed as 



(M 2 y 



2 / M 2 \ 3 

if 



(26) 



showing that for a fixed total mass M, finite-size effects will be 
smaller with a less massive secondary star. 



4.2 Mass-transfer timescale 

Next, we compare the mass-transfer timescale r^ computed from 

equations ( |l3| ) and (O) for the Roche model binaries, and equations 

(E0fl-(|24) for the numerical critical-lobe binaries. For the Roche 
— — JT\ 

model, we use Eggleton's formula in equation (|J). Note that the 

mass-radius relation ( |l0| ) is used to compute r^. In Figure Q curves 

(i) and (ii) refer to a binary system with a total mass of 0.8Mq , with 

the first one being the result of the fully numerical calculation and 
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the second one being obtained using the Roche mode(Jlj Curves (iii) 
and (iv) are similar to (i) and (ii), but refer to a system with total 
mass 1.3Mq. 

An important aspect emerging from Figure kl is that, for all of 
the relevant masses of the secondary, the mass-transfer timescale 
computed numerically is slightly larger than that obtained from the 
Roche model. In other words, the Roche model slightly underes- 
timates the evolutionary timescale, thus resulting in a larger min- 
imum orbital period. As a consequence, our results indicate that 
a treatment which refines the hydrodynamic aspects of the critical 
lobe model for CVs produces even larger discrepancies with the 
observations. 

In Figure M, we also show the curves for a parametrized 
timescale f„„ 



GMJ 



(27) 



where the adjustable coefficient a £ [0.25, 1.0] has been intro- 
duced to account for the uncertainty in the determination of the 
minimum period [cf. eq. (0)] 

For a low-mass main sequence star, the luminosity of the sec- 
ondary Z/2 a ppearing in (|27|) can be expressed as (K ippenhahn & 
Weigert 1990), 



L, 



: 0.1413 



M 2 
M 



(28) 



corresponds to the minimum 



The orbital period at which r^ = 
period of the binary system Ptum- In Tables I and II, we show the 
values of the binary parameters and the explicit values of Ptum at 
TV. = 7Vt, for several values of a and AI. 



5 CONCLUSIONS 

The Roche model has been widely used in the literature to describe 
the dynamical evolution of CVs. The minimum periods computed 
on the basis of this model have shown a disagreement with obser- 
vations and several attempts made to include more realistic descrip- 
tions of the system have not yet produced a satisfactory explanation 
for it. In order to assess the accuracy of the Roche model we have 
carried out numerical calculations of equilibria of corotating semi- 
detached binaries. The models obtained numerically are "exact" in 
the sense that they are hydrodynamically self-consistent and fully 
account for the finite size of the secondary star which fills its criti- 
cal lobe. 

However, when physical systems obtained with the two ap- 
proaches are compared, these show only minimal differences which 
are always smaller than a few percent. Therefore, our results indi- 
cate that the explanation for the mismatch between the theoretical 
minimum period and the observed one must be due to something 
else. A better agreement with observations might be found when 
mechanisms that increase the efficiency of angular momentum loss 
from the system, are fully taken into account. Among these mech- 
anisms one could consider, for instance, stellar winds or magnetic 
torques, both of which have been shown to be very efficient in re- 
moving angular momentum. An interesting alternative mechanism 
which could increase the angular momentum loss via gravitational 
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Figure 4. r^ is plotted as a function of the mass of the secondary. Curves 
(i) and (ii) are for a binary system with a total mass of 0.8A/Q. Curve 
(i) corresponds to the result of the numerical computation, whereas (ii) is 
obtained with the Roche model based on the Eggleton formula. Curves (iii) 
and (iv) are similar to (i) and (ii), but for a system with total mass 1.3Mq. 
Also plotted are the curves of f KH for the secondary. The labels of these 
curves refer to the values of the parameter a defined in equation (1271). 



model 


M 2 /M G 


q 


t m [ Gvr ] 


Ptum [min] 


a = 0.25 


Roche 
Numerical 


0.0824 
0.0821 


15.05 
15.11 


1.019 
1.025 


60.20 
60.02 


a = 0.3 


Roche 
Numerical 


0.0896 
0.0893 


13.77 
13.83 


1.1174 
1.119 


64.28 
64.08 


a = 0.35 


Roche 
Numerical 


0.0962 
0.0958 


12.77 
12.83 


1.208 
1.204 


67.96 
67.74 


a = 0.5 


Roche 
Numerical 


0.1133 
0.1127 


10.68 
10.74 


1.442 
1.449 


77.22 
76.92 


a = 1.0 


Roche 
Numerical 


0.1561 
0.1550 


7.483 
7.540 


2.034 
2.049 


99.12 
98.58 



Table 1, Comparison of binary parameters at the point along the evolution- 
ary sequence where the relation, r^ = f KH , hold. The total mass of the 
system is 1.3M0. 



I For simplicity, we do not include here the spin angular momentum con- 
tribution of the secondary in the Roche model. As shown in Figure H, this 
produces an additional error of a few percent at most. 
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Gravitational Wave Emission in CVs 



model 


M 2 /M Q 


Q 


t m [Gyr] 


Pturn [min] 


a = 0.25 


Roche 
Numerical 


0.0702 
0.0698 


10.67 
10.72 


1.215 
1.201 


53.66 
53.45 


a = 0.3 


Roche 
Numerical 


0.0763 
0.0759 


9.725 
9.785 


1.330 
1.311 


57.28 
57.05 


a = 0.35 


Roche 
Numerical 


0.0820 
0.0815 


8.997 
9.055 


1.437 
1.412 


60.58 
60.31 


a = 0.5 


Roche 
Numerical 


0.0967 
0.0960 


7.471 
7.528 


1.714 
1.727 


68.88 
68.52 


a = 1.0 


Roche 
Numerical 


0.1341 
0.1328 


5.109 
5.172 


2.401 
2.427 


88.56 
87.90 



Lasota J-R, 2000, Talk given at the International Conference "Gravitational 

Waves: A Challenge to Theoretical Astrophysics. , Trieste 
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Misner C. W., Thome K. S., Wheeler J. A., 1970, Gravitation, Freeman & 

Co., San Francisco 
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Table 2. The same as Table 1 , except that the total mass of the system is 
O.8M 



waves is given by the possible excitation of stellar oscillations in 
binary systems. These periodic perturbations, which could be ex- 
cited in binary systems such as CVs, would generate gravitational 
wave emission which would add to the one coming from the or- 
bital motion and, when r esonant, could even b e the largest source 
of gravitational radiation (Berti & Ferrari 200C ). 
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